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Mechanical  properties  and  conductivity  were  computed  for  several  composite  polymer  electrolyte 
structures.  A  multi-phase  effective  medium  approach  was  used  to  estimate  effective  conductivity.  The 
Mori-Tanaka  approach  was  applied  for  calculating  the  effective  stiffness  tensor  of  the  composites.  An 
analysis  of  effective  mechanical  properties  was  performed  in  order  to  identify  the  composite  structures, 
which  would  be  capable  of  blocking  the  dendrites  forming  in  Li-ion  battery  when  Li  metal  is  used  as  anode. 
The  stability  parameter  which  combines  both  stiffness  and  compressibility  of  the  electrolyte  was  used  in 
the  analysis.  The  calculations  were  done  over  the  wide  range  of  Young’s  modulus  of  the  polymer  matrix 
showing  the  threshold  concentration  of  the  filler  necessary  for  the  mechanical  stability.  The  results  can 
be  used  to  formulate  design  criteria  for  solid  electrolytes  that  would  exhibit  appropriate  stiffness  and 
compressibility  to  suppress  lithium  dendrite  growth  while  maintaining  high  effective  conductivities. 
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1.  Introduction 

Lithium-ion  batteries  used  in  consumer  and  automotive  appli¬ 
cations  are  traditionally  based  on  liquid  electrolytes.  Liquid 
electrolytes  are  both  thermodynamically  unstable  and  flammable, 
raising  the  concern  for  safety  of  lithium  ion  batteries.  Utilization  of 
solid  electrolytes  will  eliminate  the  need  of  use  of  flammable  sol¬ 
vents  and  will  increase  the  safety  of  Li-ion  batteries.  In  addition,  if 
solid  electrolytes  possess  appropriate  mechanical  properties,  they 
would  suppress  the  anode  roughening  and  would  prevent  Li  den¬ 
drite  growth  which  may  lead  to  shorting  the  battery.  This  property 
in  turn  will  promote  the  use  of  metallic  lithium  as  anode,  which 
will  increase  the  energy  density  of  the  battery.  It  was  recently 
postulated,  based  on  computational  results  for  the  stability  of  Li 
electrodes  [1],  that  there  is  a  threshold  for  the  shear  modulus 
above  which  the  dendritic  growth  would  be  suppressed.  The  shear 
modulus  threshold  was  found  to  decrease  with  the  compressibil¬ 
ity  of  the  electrolyte,  property  that  is  quantified  using  Poisson’s 
ratio,  v.  For  example,  for  v  =  0  (fully  compressible),  0.33,  and  0.5 
(incompressible),  a  material  with  corresponding  shear  modulus  of 
approximately  4,  2,  and  1  times  that  of  the  lithium  metal  would 
prevent  the  amplification  of  electrode  roughness. 

Invention  of  lithium  phosphorus  oxynitride  (LiPON)  solid 
electrolyte  [2-5]  for  thin  film  battery  application  opened  new 
perspectives  in  solid  state  Li-ion  technology.  LiPON  has  a  typical 
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composition  as  Li33PO3.8N0.22*  possesses  good  lithium  conductiv¬ 
ity  of  2.3  p,S  cm-1  [3,5],  and  is  stable  at  high  cell  voltage.  In  addition, 
being  a  hard  glassy  material,  LiPON  is  expected  to  suppress  lithium 
anode  roughening  and  dendrite  protrusion.  Unfortunately,  LiPON 
would  crack  under  the  tensile  stresses  developed  during  the  battery 
cycling  as  it  is  a  brittle  material. 

In  this  study,  composite  electrolyte  structures  with  conduc¬ 
tive  polymer  matrix  and  LiPON  reinforcement  are  sought  to 
overcome  this  brittleness  problem.  Typical  polymer  matrix  can 
be  represented  by  a  complex  of  lithium  salt  (LiX)  with  a  high 
molecular  weight  polymer,  typically  PEO  [6-9].  Such  combina¬ 
tion  is  promising  in  enhancing  the  overall  mechanical  properties 
while  maintaining  good  ionic  conductivity.  Significant  amount  of 
experimental  research  conducted  on  polymer  electrolytes  shows 
that  improvement  in  mechanical  properties  [10]  and  conductivity 
[6,10,11]  can  be  achieved  by  means  of  addition  of  small  particles 
to  the  polymer  matrix.  Experiments  were  conducted  mostly  with 
PEO-based  electrolytes  mixed  with  insulating  particles  of  A1203 
[6-8,11-17], Ti02  [7,15,18],  Zr02  [9],  and  Si02  [8,17].  In  all  ofthese 
cases,  a  substantial  increase  in  the  conductivity,  compared  to  the 
conductivity  of  the  polymer,  was  observed. 

It  was  suggested  that  the  main  mechanism  responsible  for 
conductivity  increase  is  the  formation  of  a  shell  consisting  of  amor¬ 
phous  polymer  around  each  particle  [  1 9-22  ].  Such  explanation  was 
based  on  the  X-ray  studies  that  revealed  a  decreased  degree  of  crys¬ 
tallinity  in  samples  with  A1203  particles  dispersed  in  a  polymer 
matrix  [23].  Similar  results  of  decrease  in  polymer  crystallinity  and 
increased  conductivity  behavior  were  obtained  from  experiments 
on  PEO  mixtures  with  higher  molecular  weight  polyacrylamide 
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Table  1 

Conductivities  (25  °C)  of  the  three  phases  of  CPE  (S  cm-1 ). 


UPON  [3,41(7/ 

(PEO)LiC104  [26]crm 

Amorphous  (PEO)LiC104  [26]crs 

2.3  xlO-6 

8.1  x  10“8 

3.59  x  10“6 

(PAAM)  [24-26].  It  has  to  be  mentioned  that  to  the  best  of  the 
author’s  knowledge  there  is  no  experimental  evidence  for  the 
highly  conductive  amorphous  polymer  layer  at  the  LiPON/Polymer 
interface.  This  is  due  to  the  fact  that  the  investigation  of  UPON  as  a 
candidate  for  composite  polymer  electrolyte  (CPE)  reinforcement 
has  just  emerged  and  the  measurement  of  properties  and/or  phase 
characterization  at  microscale  is  ongoing. 

In  the  optimized  configuration,  the  CPE  for  Li  ion  battery  with 
lithium  anode  would  be  structured  in  such  a  way  that  its  mechan¬ 
ical  properties  are  enhanced  enough  to  prevent  cracking  as  well  as 
suppress  dendrite  penetration,  while  the  conductivity  is  increased 
compared  to  the  conductivity  of  the  reinforcing  phase.  In  the 
present  investigation,  an  analysis  of  both  conductivity  and  elas¬ 
tic  properties  is  performed  for  the  composite  microstructures.  The 
reinforcement  medium  is  considered  to  be  either  spherical  par¬ 
ticles  or  flakes  of  a  given  thickness.  The  effective  conductivity  and 
effective  stiffness  tensor  are  estimated  for  different  distributions  of 
the  filler  phase.  Composites  with  conductive  reinforcement  phase 
(prototyped  after  LiPON)  are  considered.  General  effective  medium 
theory  of  McLachlan,  which  is  here  modified  for  a  three-phase 
system,  is  applied  for  determination  of  effective  composite  con¬ 
ductivity.  The  Mori-Tanaka  approach  is  used  for  estimating  the 
effective  mechanical  properties  of  the  material.  In  the  mechani¬ 
cal  analysis,  the  effect  of  the  amorphous  polymer  phase,  i.e.,  the 
highly  conductive  shell  surrounding  the  filler  components,  on  the 
mechanical  properties  of  the  CPE  is  negligible  due  to  its  very  low 
stiffness  as  compared  to  that  of  the  filler.  Finally,  several  CPE  elec¬ 
trolyte  architectures  for  Li  batteries  are  proposed  based  on  the 
results  from  the  numerical  simulations  of  CPE  properties. 

2.  Material  properties 

Material  properties  for  the  constituents  of  the  composite  elec¬ 
trolyte  were  available  from  several  references.  The  conductivities 
of  the  three  phases  at  25  °C,  which  were  obtained  from  the  corre¬ 
sponding  Arrhenius  or  VTF  equations,  are  listed  in  Table  1  together 
with  the  corresponding  references  in  which  the  equation  parame¬ 
ters  are  reported.  The  subscripts  /,  m  and  s  represent  filler,  matrix 
and  shell  phases  of  the  composite,  respectively.  Conductivity  of  the 
shell  is  assumed  to  be  equal  to  the  conductivity  of  amorphous  phase 
(PEO)  LiC104. 

Young’s  modulus  and  Poisson’s  ratio  for  LiPON  were  obtained 
from  a  recent  nano-indentation  study  [27]  where  they  are  reported 
to  be  equal  to  77  GPa  and  0.25,  respectively.  Young’s  modulus  of 
PEO  was  measured  to  be  0.2-5  GPa  for  thin  films  [28]  and  0.29  GPa 
for  highly  crystalline  bulk  material  [29].  The  Poisson’s  ratio  of  PEO 
was  0.49  [29].  The  high  value  (5  GPa)  for  the  Young’s  modulus  was 
obtained  for  thin  films  using  large  force  indentation  experiments 
[28]  and  may  be  relevant  to  the  composites  produced  by  pressing 
techniques.  For  PEO  nanofibers,  Bellan  et  al.  report  Young’s  modu¬ 
lus  values  between  0.5  and  7.0  GPa  [30].  The  chain  Young’s  modulus 
of  1 3  GPa  was  determined  from  Raman  measurements  [31]  and  can 
serve  as  the  upper  limit  of  the  elastic  modulus. 

Evaluation  criteria  for  the  mechanical  properties  of  CPE  block  is 
based  on  the  analysis  performed  by  Monroe  and  Newman  [1],  in 
which  it  was  postulated,  based  on  linear  stability  analysis,  that  Li 
surface  stability  is  achieved  when  the  shear  modulus  of  electrolyte 
would  be  greater  than  approximately  twice  the  shear  modulus  of 
lithium  metal  when  v  =  0.33.  This  shear  modulus  threshold  was 


found  to  decrease  with  Poisson’s  ratio.  A  lower  Poisson’s  ratio  - 
which  indicates  high  compressibility,  would  require  larger  shear 
modulus  for  stability  of  the  electrode  surface.  For  incompressible 
electrolytes,  i.e.,  v^0.5,  the  shear  modulus  threshold  was  esti¬ 
mated  to  be  approximately  half  as  that  for  v  =  0.33,  i.e.,  GefflGLi  =  1.0 
while  for  fully  compressible  electrolytes,  i.e.,  0,  the  shear  mod¬ 

ulus  threshold  was  approximately  twice  as  that  for  v  =  0.33,  i.e., 
Geff/Gu  =  4.0.  In  this  paper,  GLi  and  Geff  denote  the  shear  modulus  of 
lithium  metal  and  effective  shear  modulus  of  the  electrolyte  respec¬ 
tively.  Based  on  these  reported  findings  [1],  a  stability  parameter 
that  is  dependent  of  Poisson’s  ratio  is  introduced  in  this  study  by 
linear  regression  analysis  of  GeJ^/GLl  corresponding  to  the  threshold 
points  when  electrode  stability  is  achieved.  The  values  of  GefflGLi 
higher  than  the  threshold  are  considered  providing  stability  at  Li 
metal  surface  and  in  this  case  the  condition  for  stability  parameter 
(£)  is  expressed  as 

?  =  (ai  +a2v)^-  <  1  (1) 

with  fitting  constants  being:  a\  =3.503  ;  a2  =  -  5.034.  By  means  of 
Eq.  (1)  the  influence  of  both  effective  shear  modulus  and  Poisson’s 
ratio  is  considered  in  one  parameter  which  aids  further  analysis  of 
the  results. 

3.  Computational  methodology 

In  this  section,  the  computational  methodology  is  presented  for 
both  the  effective  conductivity  and  the  effective  elastic  properties 
of  the  CPE. 

3.1.  Effective  conductivity 

Schematics  of  the  particulate  CPE  microstructure,  which  is  con¬ 
sidered  for  the  effective  conductivity  calculations,  are  shown  in 
Fig.  1.  The  two-dimensional  cross  section  in  Fig.  1  can  be  represen¬ 
tative  of  a  three-dimensional  particulate  composite.  Three  distinct 
phases,  each  possessing  its  own  conductivity  and  playing  role  in 
overall  performance  of  the  electrolyte,  can  be  identified  in  Fig.  1. 
Based  on  the  conductivity  property  data  of  PEO-based  electrolytes 
and  LiPON  (Table  1 ),  the  conductivity  of  the  phases  is  ordered  in 
increasing  value  as  crf<crm<  crs.  Due  to  the  lack  of  any  experimental 
data  regarding  the  dependence  of  the  shell  thickness  on  the  parti¬ 
cle  size,  the  shell  thickness  was  considered  to  be  the  same  for  all 
the  reinforcing  elements  in  the  composite.  Since  the  existence  of 
the  highly  conductive  amorphous  polymer  layer  was  not  yet  con¬ 
firmed  at  the  LiPON/Polymer  interface,  additional  calculations  were 
performed  without  considering  the  amorphous  shell,  i.e.,  <7s  =  crm. 

Several  analytical  formulations  for  the  prediction  of  composite 
effective  conductivity  are  widely  available.  The  most  known  formu¬ 
lations  are  Brick-Layer  Model  (series  and  parallel),  Maxwell  Theory, 
McLachlan’s  Generalized  EMT  (GEMT),  and  others.  Several  repre¬ 
sentative  EMT-type  models  are  arranged  in  Table  2  [32-34].  These 
models  were  derived  for  two-phase  mixtures.  It  has  been  shown 
that  GEMT  and  Symmetric  Bruggeman  EMT  (SBMT)  models  are 
accurate  for  two-phase  particulate  composites  [32].  In  this  study, 
the  GEMT  is  extended  from  two-phase  to  three-phase  following  the 
method  suggested  by  Nan  and  Smith  [22,34].  The  series  and  paral¬ 
lel  models  are  used  to  represent  layered  structures,  which  would 
be  obtained  with  LiPON  flakes  running  parallel  or  perpendicular  to 
the  electrode  surface. 

The  two-phase  models  listed  in  Table  2  can  be  used  directly  to 
estimate  effective  conductivity  for  cases  where  the  effect  of  inter¬ 
facial  shell  layer  between  the  particle  and  the  matrix  is  neglected. 
In  order  to  take  into  account  the  different  conductivity  of  the 
shell  these  formulations  need  to  be  adjusted.  This  is  typically 
done  by  lumping  the  filler  and  shell,  i.e.,  as  another  filler  with  a 
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Fig.  1.  Schematics  of  the  3-phase  particulate  CPE:  (a)  2D  representation  of  CPE  structure  and  (b)  composite  particle  consisting  of  filler  particle  of  radius  r  and  conductive 
shell  of  thickness  8. 


conductivity  computed  from  the  conductivities  of  the  original  filler 
and  shell  [20-23].  For  instance,  in  case  of  a  spherical  particle 
(Fig.  lb),  the  effective  conductivity  (a)  of  a  composite  particle  con¬ 
sisting  of  the  filler  particle  and  conductive  shell  can  be  expressed 
following  the  MWT  approach  (d  =  3) 


a~as  =o-i)-s~af 


of  +  2oy 


Os  +  2  On 


7 


(2) 


where  the  volume  fraction  of  the  filler  in  the  composite  particle  g 
can  be  calculated  as: 


g  = 


1 


(r  +  sy  (1  +  ay 


(3) 


where  a  =  8lr  is  dimensionless  shell  thickness  parameter.  The  vol¬ 
ume  fraction  of  the  composite  particles  in  matrix  becomes  g/g 
and  the  GEMT  effective  conductivity  equation  for  the  composite 
electrolyte  is  modified  accordingly,  as: 

/g\  -(ge//)1Al 

V#/  ffl/ti  +(1  /g{f  -  lXcTejjO1^1 


(4) 


The  conductivity  of  the  three-phase  system,  which  is  shown 
in  Fig.  1,  has  two  percolation  thresholds,  g^  and  g^\  as  has 
been  shown  experimentally  [20-23].  Therefore  the  effective  con¬ 
ductivity  curve  as  a  function  of  filler  volume  fraction  has  two 
distinct  regions.  In  the  first  region,  the  conductivity  increases  due  to 
increase  in  the  fraction  of  the  highly  conductive  shell.  Equation  4  is 
valid  in  this  region.  When  the  material  is  saturated  by  the  filler  frac¬ 
tion  at  a  point  wheng  =  g  (g*  =  g),  the  microstructure  is  essentially 


Table  2 

EMT  models  for  effective  conductivity  of  a  composite  material. 


McLachlan  Generalized  Effective  Medium  Theory  (GEMT) 
Maxwell-Wagner  Theory  (MWT) 

Asymmetric  Bruggeman  Effective  Medium  Theory  (ABEM) 
Symmetric  Bruggeman  Effective  Medium  Theory  (SBEM) 
Series  Brick-Layer  Model  (SBLM) 

Parallel  Brick- Layer  Model  (PBLM) 

Landau-Lifshitz  (LL) 

Hamilton-Crosser  (HC) 


represented  by  the  shell  and  filler  phases.  Beyond  this  threshold 
point,  the  conductivity  starts  decreasing  due  to  an  increase  in  the 
volumetric  fraction  of  the  reinforcement  phase,  which  has  the  low¬ 
est  conductivity.  In  this  regime  Eq.  (4)  is  written  as 


7 


/A2-(^)Vt2 


aj^+O/g^-ncreff)^2 

(^)1/t2 


+  (1  ~g)~ 


apt2-i 


Vt2+o  M2)-ix^)1/t2 


=  0 


(5) 


The  results  presented  in  Fig.  2  indicate  that  the  proposed  three- 
phase  CPE  homogenization  model  can  be  successfully  used  to 
describe  the  experimental  data  for  effective  conductivity  for  PEO- 
Nal/Al203  [23]  and  PE0-LiC104/PAAM  [26]  composite  systems.  The 
results  shown  in  Fig.  2  were  obtained  with  effective  percolation 
thresholds  g^  and  g^2)  of  0.28  and  0.15,  respectively  [23,32]. 
The  effective  percolation  slopes  t\  and  t2  were  chosen  to  match 
the  experimental  data.  In  general,  the  higher  the  difference  in 
conductivities  of  the  matrix  and  dispersed  phase,  the  higher  the 
percolation  slope.  In  the  present  work,  it  was  chosen  q  =  1  to  mini¬ 
mize  the  amount  of  variables  involved  in  calculations  and  based  on 
the  assumption  that  the  particulate  phase  is  conductive  in  the  com¬ 
posites  considered  here.  The  results  shown  in  Fig.  2  indicate  that 
GEMT,  which  uses  experimentally  determined  percolation  thresh¬ 
olds,  is  able  to  simulate  closely  the  variation  of  the  CPE  conductivity 
as  a  function  of  the  filler  content. 


7 -r-'V'' 


<y + [  1  /ge - 1  ](creff  )T A  +[  1  /ge  - 1  ](<Tejgr  31 

aeff-af  _  M_  0m~0f 

crejj+(d--l)af  v  SJam+(d-l)af 


=  0 


g 
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af-aeff 


ejj  _i_  ( I  _  'Jejj 

af+(d- 1  )creff  '  ^  S  >  am  +(d- 1  )oeJj 


=  0 


1 


1  -g 


creff=gcrf  +  (l  —g)crm 
( creff)v 3  =g(<X/)1/3  +(  1  -g)(<xm)1/3 


ae ;j0F+(n-l)<Tm 


3  ajr+(n-l)am 


creff  -  effective  composite  conductivity:  g  -  reinforcement  volume  fraction;  ge  -  effective  percolation  threshold;  t  -  effective  percolation  slope;  d  -  system  dimensionality;  n 
-  inclusion  shape  factor. 
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Fig.  2.  Experimental  and  computed  dimensionless  effective  conductivity  using 
GEMT  model  (data  at  room  temperature). 


3.2.  Effective  elastic  constants 

Effective  stiffness  tensor  was  obtained  based  on  the 
Mori-Tanaka  approach  [35-38].  Unlike  other  homogenization 
methods,  the  Mori-Tanaka  scheme  is  applicable  to  the  composites 
with  high  volume  fractions  of  the  reinforcing  phase  since  it  intro¬ 
duces  interactions  between  the  filler  particles,  or  inclusions,  in 
the  algebraic  formulations  for  the  effective  stiffness  tensor.  In  this 
section,  this  method  is  briefly  described  and  analytical  expressions 
are  given  for  the  case  when  spherical  particles  are  embedded  in 
the  polymer  matrix. 

The  reinforcement  inclusions  in  the  composite  material  are 
treated  as  inhomogeneities.  The  main  concept  in  the  Mori-Tanaka 
approach  is  that  the  average  strain  in  inhomogeneities  equals  that 
of  a  single  inhomogeneity  embedded  in  a  matrix  subjected  to  a 
uniform  matrix  strain.  Let  denote  the  stiffness  tensor  of  the  ith 
phase  of  the  composite  material  and  L*  denote  the  effective  stiffness 
tensor  of  the  composite.  Assuming  pure  elasticity  in  each  phase 
or(i)  _  the  effective  elastic  behavior  of  the  composite  mate¬ 

rial  is  described  as  q_  =  L*e,  with  cfl\  e W  being  stress  and  strain  in 
the  phase  i  and  a,  e  being  the  overall  macroscopic  stress  and  strain 
tensors  correspondingly.  By  the  average  theorem  [35] 


UPON 


Fig.  3.  Proposed  solid  electrolyte  block  configuration  with  different  core  reinforce¬ 
ment  geometries:  (a)  particulate;  (b)  vertical  layers  and  (c)  horizontal  layers. 


Mori-Tanaka  model  and  explicitly  gives  the  result  for  the  effec¬ 
tive  stiffness  tensor  of  the  multi-phase  composite  material.  The 
effective  moduli  produced  by  Mori-Tanaka  approach  always  sat¬ 
isfy  Hashin-Shtrikman  bounds  for  two-phase  composites  [37].  An 
expression  for  the  transformation  tensor  T  can  be  found  based  on 
the  solution  to  the  problem  of  an  inclusion  embedded  into  matrix 
as: 

Ith0  =  -Lm)l  1  (9) 


where  is  the  Eshelby  tensor  derived  for  the  inclusion  i.  The 
Eshelby  tensor  is  a  function  of  the  inclusion  geometry  and  matrix 
properties.  More  detailed  information  on  the  derivation  of  Eshelby 
tensors  can  be  found  in  Desrumaux  et  al.  [36].  The  components  of 
the  Eshelby  tensor  for  spherical  inclusions  are  represented  as 


$1 122  =  $2233  =  Si  133  =  S2211  =  S33n  =  S3322  = 
Sl212  =  S2323  =  S3131  = 


Sim  =  S2222  =  S3333  = 


15(1 -u) 
7  -  5v 

15(1  -v) 


5v  —  1 
15(1  -v) 


(10) 


g  =  (6) 

the  expression  connecting  stiffness  tensors  of  the  individual  phases 
with  the  effective  stiffness  tensor  of  the  composite  material  can  be 
obtained,  as: 


Y^Si(Lw~L*)e^ 


=  0 


(7) 


If  there  is  a  transformation  law  which  connects  average  strain  in 
inhomogeneity  to  the  average  strain  of  the  host  material  (matrix) 
(e(/l))  in  form  e®  =  T(hl)e(h),  then  Eq.  (7)  can  be  transformed  into  the 
following  form 


L*  =  L(/l)  + 


Y^Si(Lw-Lw) Zm 


ghL+Y^SiliM) 


(8) 


i±h 


where  L(h)  is  the  stiffness  tensor  of  the  host  material  and 
/  is  the  identity  tensor.  Expression  in  Eq.  (8)  represents  the 


It  has  to  be  mentioned  that  in  the  present  analysis  the  effect 
of  the  amorphous  shell  surrounding  the  filler  components,  on  the 
CPE  elastic  properties  is  neglected  due  to  its  very  low  stiffness  as 
compared  to  that  of  the  filler,  effectively  resulting  in  a  two-phase 
composite  material.  The  advantage  of  the  Mori-Tanaka  scheme  is 
that  it  provides  an  explicit  formulation  for  the  effective  stiffness 
tensor. 

4.  Results  and  discussion 

Three  geometrical  configurations  that  are  typical  for  compos¬ 
ite  materials  are  evaluated  for  both  effective  conductivity  and 
mechanical  strength  using  the  methods  described  in  the  previ¬ 
ous  sections.  A  three-dimensional  design  of  the  CPE  is  considered 
for  this  purpose.  The  envisioned  sandwich  structure  for  the  CPE 
between  cathode  and  anode  is  shown  in  Fig.  3.  Two  layers  of 
LiPON  should  provide  stability  at  the  cathode  and  prevent  rough¬ 
ening  at  the  lithium  anode.  Between  them,  a  composite  structure 
is  placed,  which  consists  of  a  mixture  of  PEO-based  electrolyte, 
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Fig.  4.  Effective  particulate  CPE  conductivity  and  mechanical  stability  parameter  as 
j3  =  0.5. 

(PEO)ioLiC104,  and  UPON  reinforcement  inclusions  of  different 
geometries,  as  shown  in  Fig.  3.  From  the  computational  viewpoint, 
the  effective  properties  of  the  composite  middle  layer  are  com¬ 
puted  first  and  then  combined  with  the  two  LiPON  sheets  to  give 
the  overall  property  for  the  structure.  Reuss  and  Voight  limiting 
approximations  were  used  to  obtain  effective  properties  of  lay¬ 
ered  structures  [33].  The  effective  conductivity  is  computed  for 
y-direction  only  since  this  is  the  direction  of  interest  in  battery 
operation.  Therefore,  in  the  subsequent  sections  aeg  refers  to  effec¬ 
tive  conductivity  along  the  y-axis. 

The  results  are  organized  such  that  the  effective  conductivity 
and  mechanical  properties  can  be  shown  side-by-side  as  a  func¬ 
tion  of  the  LiPON  amount  in  the  overall  CPE.  For  the  conductivity 
results,  the  conductivities  of  the  individual  phases  were  those 
shown  in  Table  2.  Unless  otherwise  specified,  the  Young’s  mod¬ 
ulus  for  the  PEO  is  considered  to  be  5  GPa,  which  was  measured  for 
thin-films  [30],  due  to  the  thin-film  nature  of  the  solid  electrolytes. 
The  estimated  properties  are  reported  in  dimensionless  form.  The 
parameter  represents  the  ratio  of  the  thickness  of  the  two 
LiPON  sheets  over  the  total  thickness  of  the  composite  electrolyte 
block.  The  effective  conductivity  is  normalized  with  respect  to  the 


function  of  filler  volume  fraction  for  different  values  of  parameters  or.  (a)  f3  =  0.1  and  (b) 


conductivity  of  the  polymer  matrix  and  the  effective  shear  modulus 
of  the  electrolyte  block  is  normalized  with  respect  to  the  shear  mod¬ 
ulus  of  the  lithium  metal.  Effective  conductivity  is  shown  for  five 
different  values  of  the  relative  shell  thickness  a.  Cases  without  the 
conductive  shell  were  represented  by  a  =  0.  Since  it  was  assumed  in 
this  study  that  the  thickness  of  the  conductive  amorphous  shell  is 
constant,  when  present,  and  does  not  depend  on  the  particle  size, 
different  values  of  a  represent  in  fact  the  different  particle  sizes. 

The  computed  results  of  the  effective  properties  for  electrolyte 
blocks  with  a  particulate  CPE  core  (Fig.  3a)  are  represented  in  Fig.  4. 
The  results  are  reported  here  for  /I  =  0.1  (Fig.  4a)  and  0.5  (Fig.  4b). 
The  following  observations  can  be  made  based  on  the  results  in 
Fig.  4.  The  polymer  matrix  is  the  least  conductive  phase  and  the 
addition  of  the  LiPON  filler  increases  the  conductivity  of  the  CPE. 
Such  increase  obviously  happens  even  in  the  case  when  the  amor¬ 
phous  layer  is  absent  {a  =  0).  In  the  presence  of  the  conductive  shell, 
the  increase  in  the  effective  conductivity  is  followed  by  a  slight 
decrease,  corresponding  to  the  saturated  state  when  the  particles 
come  into  contact  with  each  other.  Since  the  conductivity  of  LiPON 
is  rather  close  to  the  conductivity  of  the  amorphous  PEO,  such  a 
decrease  is  not  as  dramatic  as  the  one  observed  in  mixtures  of 


Fig.  5.  Effective  conductivity  and  mechanical  stability  parameter  of  vertically  layered  composite  electrolyte  as  function  of  filler  volume  fraction  for  different  values  of 
parameters  or.  (a)  /3  =  0.1  and  (b)  f3  =  0.5. 
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non-conductive  ceramics  with  polymer  electrolytes  [23].  The  rate 
of  the  initial  increase  of  the  effective  conductivity  with  respect  to 
the  filler  volumetric  fraction  seems  to  be  proportional  to  the  value 
of  parameter  or,  i.e.,  the  maximum  conductivity  is  achieved  at  lower 
volumetric  fractions  for  larger  values  of  a.  After  the  initial  increase, 
the  curves  corresponding  to  different  values  of  parameter  a  merge 
together  and  decrease  to  the  value  characteristic  of  conductivity  of 
the  filler  phase.  Such  a  behavior  is  observed  in  both  cases  of  param¬ 
eter  p.  In  case  of  thicker  (/3  =  0.5)  UPON  boundary  sheets,  however, 
the  maximum  conductivity  achieved  is  lower  than  in  the  case  when 
the  UPON  sheet  thickness  is  0.05  times  the  total  thickness  of  the 
composite  block  (/3  =  0.1).  The  highest  enhancement  in  conductiv¬ 
ity  is  found  to  be  approximately  40  times  the  conductivity  of  PEO 
matrix.  The  arrows  in  Fig.  4  demarcate  the  region  where  the  sta¬ 
bility  parameter  §  attains  the  value  of  1.0.  It  can  be  seen  that  such 
an  enhancement  in  mechanical  properties  can  be  achieved  within 
reasonably  low  (below  random  close  packing  limit  of  0.64  [39,40]) 
volume  fractions  of  UPON :  g = 0.26  for  the  case  of  ft  =  0.1  andg  =  0.1 9 
for  the  case  of  /3  =  0.5.  Obviously,  with  thicker  UPON  sheets  sand¬ 
wiching  the  CPE,  this  criteria  would  be  met  at  the  lower  volume 
fractions  of  particulate  loading.  However,  the  difference  between 
the  threshold  volume  fractions  corresponding  to  ft  =  0.1  and  /3  =  0.5 
is  rather  small,  since,  with  the  increase  of  the  LiPON  sheet  thick¬ 
ness  comes  a  trade-off  in  terms  of  higher  Young’s  modulus  and 
lower  values  of  Poisson’s  ratio  for  the  composite  electrolyte  block, 
which  results  in  higher  value  of  §. 

Another  aspect  that  needs  to  be  investigated  is  the  gain  (or  loss) 
in  the  effective  conductivity  at  the  levels  of  the  LiPON  volume 
fractions  where  the  mechanical  stability  criterion  is  met.  When 
a  =  0.0  the  gain  in  effective  conductivity  is  slightly  higher  for  thicker 
bounding  LiPON  sheets  ( creffl<ym  =  4.1  for  /3  =  0.5  and  crejy/orm  =  3.4  for 
/3  =  0A).  In  all  other  cases  {a  j=  0.0)  the  effective  conductivity  is 
higher  when  thin  bounding  layers  of  LiPON  are  used.  This  becomes 
especially  pronounced  at  high  values  of  parameter  a  {crefflcrm  =  29 
for  /I  =  0.5  and  crefflcrm  =  39  for  /3  =  0.1  when  a  =  0.6).  Based  on  the 
above  considerations  it  can  be  stated  that  rather  thin  layers  of 
LiPON  would  provide  mechanical  properties  of  the  composite  elec¬ 
trolyte  block  that  meet  the  Li  stability  criteria  at  reasonable  volume 
fractions  of  particulate  reinforcement  without  sacrificing  the  con¬ 
ductivity  enhancement. 

The  results  for  the  vertically  layered  (VL)  CPE  (Fig.  3(b))  sand¬ 
wiched  between  the  horizontal  LiPON  sheets  are  shown  in  Fig.  5.  In 
the  case  of  layered  structures,  the  parameter  a  represents  the  ratio 


Fig.  6.  Normalized  effective  conductivity  as  a  function  of  dimensionless  shell  thick¬ 
ness  parameter  a  for  VL  composite  electrolyte. 

of  the  shell  thickness  over  the  thickness  of  the  reinforcing  layer. 
Due  to  the  directional  reinforcement  of  the  CPE  by  LiPON  layers,  the 
material  properties  cannot  be  considered  symmetric  with  respect 
to  rotation  about  y-axis  and  therefore  two  values  of  shear  modulus 
(Gxy  and  Gyz )  and  two  values  of  Poisson’s  ratio  (vxy  and  vyz)  are  con¬ 
sidered.  This  correspondingly  results  in  two  values  of  the  stability 
parameter:  §xy  and  §yz.  The  electrode  stability  criteria  are  applied 
to  both  values.  It  can  be  seen  from  Fig.  5  that  while  the  §yz  =  1.0 
criterion  is  met  at  low  volume  fraction  of  LiPON  layers,  the  same 
criteria  for  §xy  is  satisfied  at  rather  high  content  of  reinforcement 
(~80%  and  ~60%  for  ft  =  0.1  and  0.5  respectively).  It  should  be  men¬ 
tioned  that  such  high  volume  fractions  are  practically  achievable  in 
case  of  layered  composites.  Obviously,  the  stability  parameter  cor¬ 
responding  to  the  stiffer  direction  yz  does  not  change  significantly 
with  increase  of  the  thickness  of  boundary  LiPON  layers. 

Behavior  of  effective  conductivity  at  reinforcement  volume  frac¬ 
tions  when  the  §  =  1.0  stability  criterion  is  met  is  shown  in  Fig.  6. 
The  normalized  effective  conductivity  is  shown  as  a  function  of  con¬ 
ductive  shell  thickness  parameter  a  for  two  values  of  parameter  p. 
This  way  the  effect  of  increasing  the  thickness  of  boundary  LiPON 


Fig.  7.  Effective  conductivity  and  mechanical  properties  of  horizontally  layered  composite  electrolyte  as  function  of  filler  volume  fraction  for  different  values  of  parameters 
a:  (a)  =  0.1  and  (b)  =  0.5. 
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layers  on  effective  conductivity  can  be  estimated.  Because  the 
§xy  =  1.0  condition  is  met  at  higher  volume  fractions,  the  effective 
conductivity  does  not  change  significantly  with  increase  in  both 
a  and  At  the  same  time  the  normalized  effective  conductivity 
increases  by  about  factor  of  1.5  at  high  values  of  a  when  thicker 
boundary  sheets  are  used  (Fig.  6).  It  should  be  mentioned  that  in 
order  to  suppress  amplification  of  electrode  surface  roughness,  the 
stability  criterion  should  be  met  in  both  directions  for  such  trans¬ 
versely  isotropic  material  and  therefore  the  composite  is  expected 
to  be  fabricated  with  high  volume  concentration  of  UPON  layers. 
The  results  for  the  horizontally  layered  (HL)  CPE  (Fig.  3(c))  are 
shown  in  Fig.  7.  The  mechanical  properties  are  symmetric  with 
respect  to  rotation  about  y-axis  and  therefore  Gxy  =  Gyz  =  Geff  and 
vxy  =  Vyz  =  veff.  One  stability  criterion  §  =  1.0  thus  is  used.  It  can  be 
seen  from  Fig.  7  that  high  content  of  LiPON  is  needed  to  main¬ 
tain  stability  of  the  electrode  surface  when  horizontally  layered 
composite  electrolyte  is  used.  However,  such  high  concentrations 
are  technically  achievable  and  provide  rather  good  enhancement  in 
effective  conductivity  when  amorphous  conductive  shell  is  absent 
(a  =  0). 

Analysis  of  the  results  shown  in  Figs.  4-7  for  three  different 
large-scale  reinforcement  geometries  reveals  the  interplay  of  dif¬ 
ferent  effective  properties  that  make  a  composite  electrolyte  block 
suitable  for  battery  applications  when  the  lithium  metal  anode  is 
employed.  The  configuration  suggested  here,  with  the  placement 
of  CPE  between  two  LiPON  layers,  was  sought  to  improve  the  stiff¬ 
ness  of  electrolyte  and  further  prevent  the  anode  roughening  and 
dendrite  formation.  For  some  configurations,  the  analysis  is  more 
complicated  as  more  variables  had  to  be  considered  and  two  sta¬ 
bility  criteria  must  be  satisfied.  In  some  cases,  the  sandwich  CPE 
may  reduce  the  benefit  of  improved  conductivity  of  the  electrolyte. 
The  higher  values  of  shear  modulus  are  achieved  by  either  increas¬ 
ing  the  thickness  of  boundary  LiPON  layers  or  by  increasing  the 
reinforcement  content  in  the  middle  CPE  section.  In  both  cases  this 
results  in  lower  values  of  Poisson’s  ratio.  The  combined  influence 
of  the  shear  modulus  and  Poisson’s  ratio  on  the  stability  of  the  elec¬ 
trode  surface  can  be  assessed  by  the  stability  parameter  §  (Eq.  (1 )) 
which  was  deduced  from  the  reported  results  of  effect  of  mechani¬ 
cal  properties  of  electrolyte  on  amplification  of  electrode  roughness 
[1]. 

Figs.  4-7  demonstrate  that  the  stability  criteria  can  be  satis¬ 
fied  in  all  of  the  three  cases  of  reinforcement  geometry  inside 
the  CPE  layer.  The  lowest  reinforcement  volume  fraction,  which 


Fig.  8.  Effect  of  the  polymer  Young’s  modulus  ( EP )  on  the  stability  parameter  for 
particulate  CPE. 

allows  for  §  =  1.0,  is  found  in  CPE  with  particulate  reinforcement 
and  LiPON  sheet  thickness  of  25%  of  the  total  thickness  of  the  elec¬ 
trolyte  (Fig.  4b).  It  should  be  re-iterated  here  that  the  existence 
of  amorphous  polymer  interfacial  layer  has  not  been  experimen¬ 
tally  confirmed  for  LiPON/PEO  systems.  Therefore  the  results  for 
a  ±  0  are  presented  for  comparison  purposes  and  the  discussion 
should  be  mainly  based  on  the  behavior  of  mixture  of  bare  LiPON 
with  PEO  electrolyte.  From  this  perspective,  the  highest  conductiv¬ 
ity  at  point  when  §  =  1 .0  is  achieved  in  vertically  layered  electrolyte 
blocks,  with  the  conductivity  being  improved  by  22  times  as  com¬ 
pared  to  polymer  matrix  conductivity  (Fig.  6).  This  improvement  is 
achieved  at  rather  large  content  of  LiPON  layers  (~80%).  However 
the  benefit  of  conductivity  improvement  could  make  this  arrange¬ 
ment  the  most  promising.  The  20%  of  the  volume  in  this  case  would 
be  occupied  by  the  polymer,  which  may  improve  the  ductility  of  the 
electrolyte  and  prevent  it  from  cracking. 

For  the  sake  of  completeness,  results  for  the  stability  thresh¬ 
old  parameter  are  now  presented  for  other  values  of  the  Young’s 
modulus  of  the  polymer.  Fig.  8  shows  the  effect  of  polymer  elastic 
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Fig.  9.  Amount  of  reinforcement  required  to  achieve  mechanical  stability  of  CPE  as  a  function  of  polymer  Young’s  modulus:  (a)  particulate  CPE  and  (b)  VL  and  HL  CPE. 
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modulus  (Ep)  on  the  behavior  of  stability  parameter  for  the  case  of 
particulate  CPE.  The  values  of  the  modulus  range  from  0.29  GPa  to 
the  chain  limit  of  13  GPa  determined  from  Raman  measurements 
[31].  The  points  where  stability  criterion  is  satisfied  are  marked 
with  asterisks.  Apparently  rather  high  loadings  of  filler  are  required 
to  achieve  stability  as  the  modulus  approaches  its  lower  limit. 

The  volume  fractions  (g*)  corresponding  to  §  =  1.0  are  plotted 
against  the  Young’s  modulus  of  the  polymer  in  Fig.  9.  Since  the  HL 
composite  represents  the  90°  rotation  of  the  VL  composite  about 
z-axis,  which  is  the  axis  of  symmetry  for  transversely  isotropic 
material,  the  values  of  §  in  xy  plane  coincide  for  both  VL  and  HL  cases 
and  are  represented  in  one  plot  (Fig.  9b).  Interestingly,  the  effect  of 
increase  in  boundary  reinforcement  sheets  is  not  pronounced  for 
the  case  of  particulate  CPE  (Fig.  9a).  It  can  be  seen  from  Fig.  9  that  no 
reinforcement  is  necessary  to  achieve  the  stability  when  the  poly¬ 
mer  elastic  modulus  reaches  values  of  9  GPa  and  1 1  GPa  for  /3  =  0.5 
and  0.1  respectively.  These  values  are  rather  high  even  for  the 
most  stiff  block  co-polymers  synthesized  nowadays  for  Li  ion  bat¬ 
tery  electrolyte  [41,42].  Therefore,  it  can  be  stated  that  fabrication 
of  composite  systems  is  necessary  in  order  to  suppress  electrode 
roughening.  Particle  reinforced  isotropic  composites  seem  to  be  the 
best  choice  to  achieve  the  stability  with  lower  amount  of  stiffening 
phase.  The  results  reported  in  this  study  can  be  used  as  guidelines 
for  CPE  design.  Finally,  the  choice  of  geometrical  configurations  for 
CPE  electrolytes  would  also  depend  on  the  manufacturing  require¬ 
ments  to  produce  such  an  arrangement  with  minimal  cost. 

5.  Conclusions 

The  data  presented  in  this  study  for  three  different  large-scale 
reinforcement  geometries  reveals  the  interplay  among  several 
effective  properties  that  make  a  composite  electrolyte  block  suit¬ 
able  for  battery  applications  when  the  lithium  metal  anode  is 
employed.  A  homogenization  model  for  a  three-phase  media  was 
used  to  evaluate  the  effective  CPE  conductivity.  Effective  stiff¬ 
ness  tensor,  from  which  the  shear  modulus  and  Poisson’s  ratio 
were  estimated,  was  obtained  based  on  the  Mori-Tanaka  approach. 
The  dependence  of  effective  conductivity  and  electrode  mechani¬ 
cal  stability  parameter  on  such  variables  as  filler  volume  fraction, 
thickness  of  highly  conductive  interfacial  layer,  and  polymer  elastic 
modulus  has  been  studied.  Both  elastic  and  compressible  proper¬ 
ties,  through  shear  modulus  and  Poisson’s  ratio,  are  considered 
to  formulate  a  mechanical  stability  criterion.  In  order  to  increase 
the  resistance  to  electrode  roughening,  the  proposed  electrolyte 
structure  consisted  of  CPE  (either  particulate  or  layered)  block 
sandwiched  between  two  horizontal  layers  of  LiPON.  The  results 
indicate  that  the  lowest  volume  fraction  of  reinforcement  capa¬ 
ble  of  making  electrolyte  suitable  for  dendrite  suppression  can  be 
achieved  with  particulate  CPE.  At  the  same  time  the  best  conduc¬ 
tivity  would  be  achieved  in  vertically  layered  structures  at  high 
reinforcement  content. 
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